• 专利附录
  • 专利说明
  • 权利要求
  • 类似技术
  • 同族专利
  • 引用文献
  • 法律状态
  • 优先权
    • 专利摘要:
    The present disclosure discloses a magnetic compensation method based on an aeromagnetic compensation error model, comprising: acquiring an upper limit of an error of a magnetic noise by using the aeromagnetic compensation error model, before an aeromagnetic flight, wherein the magnetic noise is caused by both an environmental background field in an exploration area and an aeromagnetic flight platform; determining, according to the upper limit of the error of the magnetic noise, whether the environmental background field and the aeromagnetic flight platform are suitable for an aeromagnetic survey flight, and if so, performing a calibration flight to acquire a compensation coefficient; and acquiring data of an attitude term by performing an actual flight, calculating an interference magnetic field by the data of the attitude term and the compensation coefficient acquired during the calibration flight, and performing magnetic compensation. This method may may estimate an environmental background field in an exploration area before an aeromagnetic compensation calibration flight is performed, and magnetic characteristics of an aeromagnetic exploration platform are estimated effectively in advance, to determine in advance whether the exploration area and an aeromagnetic flight platform are suitable for aeromagnetic exploration, which provides a theoretical basis for the selection and construction of the platform, and avoids occurrence of an invalid flight.
    公开号:EP3690491A1
    申请号:EP17926970.9
    申请日:2017-09-30
    公开日:2020-08-05
    发明作者:Qunying ZHANG；Peilin Wu；Luzhao CHEN；Guangyou FANG；Wanhua ZHU
    申请人:Institute of Electronics of CAS；
    IPC主号:G01V3-00
    专利说明:
    [0001] The present disclosure belongs to a field of geophysical aeromagnetic exploration, and relates to a magnetic compensation method based on an aeromagnetic compensation error model. BACKGROUND
    [0002] Aeromagnetic exploration, as an important aerial geophysical survey method, has been widely used in the field of geophysics. The exploration method has significant advantages such as high efficiency, safety, reliability etc. In an aeromagnetic exploration process, a device carrying platform is usually a fixed wing, a helicopter, etc., in which an optically pumped magnetometer is usually fixed on the carrying platform. Because an aircraft is made of a ferromagnetic material, in a geomagnetic field environment, the aircraft itself may generate an interference magnetic field signal to a signal detected by the optically pumped magnetometer. This interference magnetic field may cause failure to extract valid magnetic anomaly information from magnetic measurement data. Therefore, it is of great significance to effectively compensate for the interference magnetic field of the aircraft.
    [0003] At present, foreign magnetic compensation devices, for example, AADC series magnetic compensators from the RMS company and aeromagnetic compensation devices from the PICO company, are mainly used in the field of aeromagnetic exploration in domestic. A compensation algorithm for the above compensation devices is based on the conventional aeromagnetic compensation algorithm design, and is characterized in that data collected by an optically pumped magnetometer and a fluxgate magnetometer is firstly passed through a low-pass filter to remove some noises which are not related to a magnetic interference of the aircraft therefrom, and then is subjected to the least square algorithm to achieve the purpose of removing an interference magnetic field of the aircraft.
    [0004] However, the conventional aeromagnetic compensation model may not be used to determine an upper limit of compensation of its own, and may not predict in advance whether an environmental background field in an area to be explored satisfies requirements of aeromagnetic exploration. In addition, in the existing aeromagnetic compensation models, there is no way to select a theoretical determination method suitable for a flight platform, and it may not be predicted whether magnetic characteristics of the flight platform satisfy the requirements of aeromagnetic exploration before an aeromagnetic calibration flight is performed, and therefore, it may cause an invalid flight to be performed in an actual aeromagnetic operation. SUMMARY(1) Technical problem to be solved
    [0005] In view of the above problem, the present disclosure provides a magnetic compensation method based on an aeromagnetic compensation error model to predict whether an area to be explored, an environmental background field and a flight platform are suitable for aeromagnetic exploration, thereby avoiding an invalid flight. (2) Technical solutions
    [0006] According to an aspect of the present disclosure, the present disclosure discloses a magnetic compensation method based on an aeromagnetic compensation error model, comprising: acquiring an upper limit of an error of a magnetic noise by using the aeromagnetic compensation error model, before an aeromagnetic flight, wherein the magnetic noise is caused by both an environmental background field in an exploration area and an aeromagnetic flight platform; determining, according to the upper limit of the error of the magnetic noise, whether the environmental background field and the aeromagnetic flight platform are suitable for an aeromagnetic survey flight, and if so, performing a calibration flight to acquire a compensation coefficient; and acquiring data of an attitude term by performing an actual flight, calculating an interference magnetic field by the data of the attitude term and the compensation coefficient acquired during the calibration flight, and performing magnetic compensation.
    [0007] In some embodiments of the present disclosure, the aeromagnetic compensation error model is established by: acquiring various interference magnetic fields, and acquiring errors of magnetic interferences corresponding to the interference magnetic fields based on expressions of the respective interference magnetic fields, so as to acquire a total magnetic interference error.
    [0008] In some embodiments of the present disclosure, the various interference magnetic fields comprise: a constant interference magnetic field, an induced interference magnetic field, and an eddy current interference magnetic field, and the total magnetic interference error in the aeromagnetic compensation error model satisfies:Error T = error P + error I + error E
    [0009] HPERM (t) is the constant interference magnetic field; HIND (t) is the induced interference magnetic field; HEDDY (t) is the eddy current interference magnetic field; T(t), L(t) and V(t) are outputs of a three-axis fluxgate magnetometer in three directions; and ΔT, ΔV and ΔL are error levels of measured values of the fluxgate in three directions respectively.
    [0010] In some embodiments of the present disclosure, the constant interference magnetic field HPERM (t) satisfies:H PERM t = c 1 u 1 + c 2 u 2 + c 3 u 3= A 1 c 1 + A 2 c 2 + A 3 c 3 ,
    [0011] In some embodiments of the present disclosure, the direction cosines u 1, u 2 and u 3 satisfy:u 1 = cos X t = T t ⁄ T t 2 + L t 2 + V t 2
    [0012] In some embodiments of the present disclosure, determining, according to the upper limit of the error of the magnetic noise, whether the environmental background field and the aeromagnetic flight platform are suitable for an aeromagnetic survey flight comprises: determining whether the environmental background field and the aeromagnetic flight platform satisfy an aeromagnetic experiment, according to whether a calculated aeromagnetic compensation upper limit Is less than a lower limit of accuracy for the aeromagnetic experiment, and if the aeromagnetic compensation upper limit is less than the lower limit of the accuracy for the aeromagnetic experiment, determining that it is suitable for the aeromagnetic survey flight; otherwise, determining that it is not suitable for the aeromagnetic survey flight.
    [0013] In some embodiments of the present disclosure, performing a calibration flight to acquire a compensation coefficient comprises: performing the calibration flight and acquiring an attitude term according to outputs of a three-axis fluxgate magnetometer in three directions and strength of a geomagnetic field; acquiring an interference magnetic field of an aircraft at an optically pumped magnetometer according to data of the calibration flight; and calculating the compensation coefficient using the attitude term and the interference magnetic field, whereinthe compensation coefficient satisfies: C =A T A − 1A T H d
    [0014] In some embodiments of the present disclosure, the interference magnetic field of the aircraft at the optically pumped magnetometer satisfies:H d t = H PERM t + H IND t + H EDDY t
    [0015] In some embodiments of the present disclosure, a magnetic field after the magnetic compensation is performed satisfies:H C = H M − A M C
    [0016] In some embodiments of the present disclosure, acquiring the upper limit of the error of the magnetic noise caused by the environmental background field in the exploration area comprises: measuring the environmental background magnetic field at a certain point in the area to be explored in a state without an exploration aircraft, and acquiring the upper limit of the error of the magnetic noise caused by the environmental background field, by using the aeromagnetic compensation error model; and acquiring the upper limit of the error of the magnetic noise caused by the aeromagnetic flight platform comprises: measuring a magnetic field of the noise in an operation state, wherein, in the operation state, devices such as engines and avionics of an aircraft are turned on during an exploration operation, and acquiring the upper limit of the error of the magnetic noise caused by the aeromagnetic flight platform, by using the aeromagnetic compensation error model. (3) Beneficial effects
    [0017] As may be seen from the above technical solutions, the magnetic compensation method based on an aeromagnetic compensation error model according to the present disclosure has the following beneficial effects.
    [0018] Before an aeromagnetic compensation flight is performed, a model error caused by an environmental background magnetic noise, engines of an aircraft and avionics is firstly calculated to determine an upper limit of a background error, so as to determine whether a magnetic noise level in an area or of the aircraft is suitable for an aeromagnetic survey flight. There is provided a method for calculating an upper limit of magnetic compensation, which may estimate an environmental background field in an exploration area before an aeromagnetic compensation calibration flight is performed, and magnetic characteristics of an aeromagnetic exploration platform are estimated effectively in advance, to determine in advance whether the exploration area and an aeromagnetic flight platform are suitable for an aeromagnetic exploration, which provides a theoretical basis for the selection and construction of the platform, and avoids occurrence of an invalid flight. BRIEF DESCRIPTION OF THE DRAWINGS
    [0019] Fig. 1 is a flowchart of a magnetic compensation method based on an aeromagnetic compensation error model according to an embodiment of the present disclosure. Fig. 2 is a calculation block diagram of an aeromagnetic compensation error model according to an embodiment of the present disclosure. Fig. 3 is a diagram of output results of three channels of a fluxgate under a background field in a measurement area according to an embodiment of the present disclosure. Fig. 4 is error curves of three types of magnetic interference and a total magnetic interference under a background field in a measurement area according to an embodiment of the present disclosure. Fig. 5 is a diagram of output results of three channels of a fluxgate in an experiment using an automobile to simulate magnetic noises of engines of an aircraft and avionics according to an embodiment of the present disclosure. Fig. 6 is error curves of three types of magnetic interference and a total magnetic interference in an experiment using an automobile to simulate magnetic noises of engines of an aircraft and avionics according to an embodiment of the present disclosure.DETAILED DESCRIPTION
    [0020] The present disclosure provides a magnetic compensation method based on an aeromagnetic compensation error model. Before an aeromagnetic compensation flight is performed, a model error caused by an environmental background magnetic noise, engines of an aircraft and avionics is firstly calculated to determine an upper limit of a background error, so as to determine whether a magnetic noise level in an area or of the aircraft is suitable for an aeromagnetic survey flight. There is provided a method for calculating an upper limit of magnetic compensation, which may estimate an environmental background field in an exploration area before an aeromagnetic compensation calibration flight is performed, and magnetic characteristics of an aeromagnetic exploration platform are estimated effectively in advance, to determine in advance whether the exploration area and an aeromagnetic flight platform are suitable for an aeromagnetic exploration, which provides a theoretical basis for the selection and construction of the platform, and avoids occurrence of an invalid flight.
    [0021] In order to make the purposes, technical solutions, and advantages of the present disclosure more clear, the present disclosure will be further described in detail below in combination with specific embodiments and with reference to accompanying drawings.
    [0022] In a first exemplary embodiment of the present disclosure, there is provided a magnetic compensation method based on an aeromagnetic compensation error model.
    [0023] Fig. 1 is a flowchart of a magnetic compensation method based on an aeromagnetic compensation error model according to an embodiment of the present disclosure. Fig. 2 is a calculation block diagram of an aeromagnetic compensation error model according to an embodiment of the present disclosure.
    [0024] As shown in Figs. 1 and 2, the magnetic compensation method based on an aeromagnetic compensation error model according to the present disclosure comprises the following steps.
    [0025] In step S102, the aeromagnetic compensation error model is established.
    [0026] As shown in Fig. 2, in this step, establishing an aeromagnetic compensation error model comprises: acquiring various interference magnetic fields, and acquiring errors of magnetic interferences corresponding to the interference magnetic fields based on expressions of the respective interference magnetic fields, to acquire a total magnetic interference error, so as to establish the aeromagnetic compensation error model.
    [0027] In the conventional aeromagnetic exploration process, various interference magnetic fields are mainly divided into three types of interference magnetic fields, comprising: a constant interference magnetic field, an induced interference magnetic field and an eddy current interference magnetic field. Therefore, for the three types of interference magnetic fields, magnetic interference errors corresponding to the interference magnetic fields are derived and comprise: an error of a constant magnetic interference, an error of an induced magnetic interference and an error of an eddy current magnetic interference, so as to obtain a total magnetic interference error.
    [0028] Expressions of the constant interference magnetic field, the induced interference magnetic field and the eddy current interference magnetic field in aeromagnetic flight data will be described below, and the corresponding errors and the total error are derived accordingly.
    [0029] Firstly, an attitude of an aircraft is measured by establishing a relationship between outputs of a three-axis fluxgate magnetometer and direction cosines of a geomagnetic field vector in an aircraft carrier coordinate system. The relationship therebetween may be expressed as:u 1 = cos X t = T t ⁄ T t 2 + L t 2 + V t 2 u 2 = cos Y t = L t ⁄ T t 2 + L t 2 + V t 2 u 3 = cos Z t = V t ⁄ T t 2 + L t 2 + V t 2
    [0030] An interference magnetic field generated by the aircraft at an optically pumped magnetometer is composed of the constant interference magnetic field, the induced interference magnetic field, and the eddy current interference magnetic field. The interference magnetic field of the aircraft at the optically pumped magnetometer is expressed as:H d t = H PERM t + H IND t + H EDDY t= ∑ i = 1 18A i c i
    [0031] A corresponding matrix expression is expressed as:H d = A C + z
    [0032] In this implementation, the observation noise is a Gaussian white noise with a mean of 0.wherein the constant interference magnetic field is expressed as follows:H PERM t = c 1 u 1 + c 2 u 2 + c 3 u 3= A 1 c 1 + A 2 c 2 + A 3 c 3
    [0033] The induced interference magnetic field is expressed as follows:H IND t = H e tc 4 u 1 2 + c 5 u 1 u 2 + c 6 u 1 u 3+ c 7 u 2 2 + c 8 u 2 u 3 + c 9 u 3 2 = A 4 c 4 + A 5 c 5 + A 6 c 6 + A 7 c 7 + A 8 c 8 + A 9 c 9
    [0034] The eddy current interference magnetic field is expressed as follows:H EDDY t = H ec 10 u 1 u 1 ′ + c 11 u 1 u 2 ′ + c 12 u 1 u 3 ′ + c 13 u 2 u 1 ′ + c 14 u 2 u 2 ′+ c 15 u 2 u 3 ′ + c 16 u 3 u 1 ′ + c 17 u 3 u 2 ′ + c 18 u 3 u 3 ′ = A 10 c 10 + A 11 c 11 + A 12 c 12 + A 13 c 13 + A 14 c 14 + A 15 c 15 + A 16 c 16 + A 17 c 17 + A 18 c 18
    [0035] Secondly, the corresponding errors and the total error are derived based on the expressions of the constant interference magnetic field, the induced interference magnetic field and the eddy current interference magnetic field.
    [0036] The error of the constant magnetic interference may be obtained according to the expression of the constant interference magnetic field as shown in formula (8) in combination with a measurement relationship of the attitude of the aircraft as shown in formula (1), and is expressed as follows:error P = ∂ H PERM∂ T Δ T + ∂ H PERM∂ L Δ L + ∂ H PERM∂ V Δ V
    [0037] Partial derivatives of the constant interference magnetic field with respect to the measured values in the three directions satisfy the following relationship: ∂ H PERM∂ T = c 1 L 2 + V 2 H e 3 − c 2 LT H e 3 − c 3 VT H e 3
    [0038] The error of the induced magnetic interference may be obtained according to the expression of the induced interference magnetic field as shown in formula (9) in combination with the measurement relationship of the attitude of the aircraft as shown in formula (1), and is expressed as follows:error I = ∂ H IND∂ T Δ T + ∂ H IND∂ L Δ L + ∂ H IND∂ V Δ V
    [0039] Partial derivatives of the induced interference magnetic field with respect to the measured values in the three directions satisfy the following relationship:∂ H IND∂ T = c 4 T 3 + 2 TL 2 + 2 TV 2 H e 3 + c 5 L 3 + V 2 LH e 3 + c 6 V 3 + L 2 VH e 3 − c 7TL 2H e 3 − c 8 LVT H e 3 − c 9 V 2 TH e 3
    [0040] The error of the eddy current magnetic interference may be obtained according to the expression of the eddy current interference magnetic field as shown in formula (10) in combination with the measurement relationship of the attitude of the aircraft as shown in formula (1), and is expressed as follows:error E = ∂ H EDDY∂ T Δ T + ∂ H EDDY∂ L Δ L + ∂ H EDDY∂ V Δ V
    [0041] Partial derivatives of the eddy current interference magnetic field with respect to the measured values in the three directions satisfy the following relationship:∂ H EDDY∂ T = c 10 T H e′ + T L 2 + V 2H e 3′ + c 11 L H e′ − TTL H e 3′ + c 12 V H e′ − TTV H e 3′ + c 13 L L 2 + V 2H e 3′ − c 14 LTL H e 3′ − c 15 LTV H e 3′+ c 16 V L 2 + V 2H e 3′ − c 17 VTL H e 3′ − c 18 VTV H e 3′
    [0042] According to the obtained error of the constant magnetic interference, the obtained error of the induced magnetic interference and the obtained error of the eddy current magnetic interference, the total magnetic interference error is obtained, which is the aeromagnetic compensation error model. Here, the total magnetic interference error satisfies:Error T = error P + error I + error E
    [0043] ΔT, ΔV and ΔL are obtained using the measured values of the attitude of the aircraft, and the aeromagnetic compensation error model may be obtained according to the formula (23).
    [0044] In step S104, an upper limit of an error of a magnetic noise is acquired by using the aeromagnetic compensation error model, before an aeromagnetic flight, wherein the magnetic noise is caused by both an environmental background field in an exploration area and an aeromagnetic flight platform.
    [0045] Firstly, the error caused by the magnetic noise existing in the environmental background field in the exploration area is acquired. In a state where there is no exploration aircraft, the environmental background magnetic field is measured at a certain point in the area to be explored, and the error of the magnetic noise caused by the environmental background field is acquired by using the aeromagnetic compensation error model.
    [0046] In this embodiment, an experimental scene is as follows: an open field in the wild is selected as a measurement area, an optically pumped magnetometer and a fluxgate magnetometer are fixed on an extension rod made of a nonmagnetic material, the extension rod is fixed to a wooden stand bracket which is located at a certain point in the measurement area, and it is ensured that there is no uncontrollable magnetic noise source such as motor automobiles in a surrounding environment, to ensure that the measured data is a magnetic noise only caused by an environmental background field.
    [0047] A process of processing experimental data comprises: processing test results of the fluxgate magnetometer to obtain error levels of the measured values of the fluxgate in three directions, and obtaining errors of the magnetic noise caused by the environmental background field based on expressions of respective components in the aeromagnetic compensation error model.
    [0048] Fig. 3 is a diagram of output results of three channels of a fluxgate under a background field in a measurement area according to an embodiment of the present disclosure. Fig. 4 is error curves of three types of magnetic interference and a total magnetic interference under a background field in a measurement area according to an embodiment of the present disclosure.
    [0049] In the present embodiment, data of a fluxgate magnetometer in three directions is processed by a high-pass filter of 0.065Hz to obtain errors ΔT, ΔV and ΔL of measured values of the fluxgate in the three directions, and results of ΔT, ΔV and ΔL from which direct current components are removed are shown in (a), (b) and (c) in Fig. 3 respectively. Error curves of three types of magnetic interference and a total magnetic interference under a background field in a measurement area, shown in Fig. 4, are obtained according to the formulas (23), (11), (15) and (19), wherein (a) is the error curve of the constant magnetic interference, (b) is the error curve of the induced magnetic interference, (c) is the error curve of the eddy current magnetic interference, and (d) is the error curve of the total magnetic interference.
    [0050] Table 1 shows standard deviations and Figure Of Merit (FOM) values of the three types of magnetic interference fields and the total magnetic interference field under the environmental background field, and standard deviations and FOMs of different magnetic interference fields are calculated from the curves shown in Fig. 4 and are shown in Table 1.
    [0051] Table 1 Standard deviations and FOMs of different magnetic interference fields in environmental background field Type Standard deviation(nT) FOM(nT) Constant magnetic interference 2.85×10-7 3.04×10-6Induced magnetic interference 8.80×10-3 7.45×10-2Eddy current magnetic interference 3.73×10-6 2.78×10-4Total magnetic interference field 8.80×10-3 7.45×10-2
    [0052] It may be known from Table 1 that when an aeromagnetic exploration experiment is carried out in this measurement area, no matter what kind of aircraft is used as a magnetic measurement platform, an upper limit of aeromagnetic compensation which may be achieved is 7.45×10-2nT. Thus, this area has a good background magnetic environment, and is suitable for a high-accuracy aeromagnetic operation.
    [0053] Secondly, the error caused by the magnetic noise existing in an aeromagnetic flight platform is acquired. In the aeromagnetic flight platform, main sources of errors in magnetic noises are engines of an aircraft and avionics. Therefore, during an exploration operation, the engines of the aircraft and the avionics are in a turn-on state, a noise magnetic field in this operation state is measured, and an error in a magnetic noise caused by the aeromagnetic flight platform is acquired by using the aeromagnetic compensation error model.
    [0054] Since the magnetic noises of the engines of the aircraft and the avionics are generated in a similar principle to that of magnetic noises of engines of automobiles and on-board electronic devices to a certain extent, in the present embodiment, an automobile is used to simulate the magnetic noise generated by the aircraft and measurement is made by way of experimental simulation in the wild.
    [0055] In this embodiment, an experimental scene is as follows: an open field in the wild is selected as a measurement area, an optically pumped magnetometer and a fluxgate magnetometer are fixed on an extension rod made of a nonmagnetic material, the extension rod is fixed to a wooden stand bracket which is located at a certain point in the measurement area, a modern SUV is used to simulate an aircraft platform and is at an idle speed during the experiment, and on-board electronic devices are in a turn-on state.
    [0056] Magnetic noises in the experimental data obtained in this experimental scene comprise an error of a magnetic noise caused by an environmental background field and an error of magnetic noises caused by engines of the automobile and the on-board electronic devices.
    [0057] Fig. 5 is a diagram of output results of three channels of a fluxgate in an experiment using an automobile to simulate magnetic noises of engines of an aircraft and avionics according to an embodiment of the present disclosure. Fig. 6 is error curves of three types of magnetic interference and a total magnetic interference in an experiment using an automobile to simulate magnetic noises of engines of an aircraft and avionics according to an embodiment of the present disclosure.
    [0058] Data of a fluxgate magnetometer in three directions is processed according to a method of acquiring a magnetic noise existing in an environmental background field to obtain errors ΔT, ΔV and ΔL of measured values of the fluxgate in the three directions, and results of ΔT, ΔV and ΔL from which direct current components are removed are shown in (a), (b) and (c) in Fig. 5 respectively. Error curves of three types of magnetic interference and a total magnetic interference caused by an aeromagnetic flight platform and a background field in an observation area during an exploration process are obtained according to the formulas (23), (11), (15) and (19) and are shown in Fig. 6, wherein (a) is the error curve of the constant magnetic interference, (b) is the error curve of the induced magnetic interference, (c) is the error curve of the eddy current magnetic interference, and (d) is the error curve of the total magnetic interference.
    [0059] It may be seen from Fig. 5 that there is a step-like fluctuation in the measured signal. The fluctuation signal is mainly caused by automatic gear shifting of a gearbox of a SUV, and belongs to a sudden signal.
    [0060] Table 2 shows standard deviations and FOM values of the three types of magnetic interference fields and the total magnetic interference field in an exploration mode, and standard deviations and FOMs of different magnetic interference fields are calculated from the curves shown in Fig. 6 and are shown in Table 2.
    [0061] Table 2 Standard deviations and FOMs of different magnetic interference fields in the exploration mode Type Standard deviation (nT) FOM(nT) Constant magnetic interference 5.66×10-6 9.92×10-5Induced magnetic interference 0.49 6.76 Eddy current magnetic interference 2.83×10-4 1.26×10-2Total magnetic interference field 0.49 6.76
    [0062] It may be seen by comparing Table 2 with Table 1 that in the environmental background field, an theoretical upper limit of compensation which may be achieved by the aeromagnetic compensation model is ideal, and is only 7.45×10-2 nT. In the exploration mode, a total upper limit is 6.76 nT, which indicates that the magnetic noises caused by the engines of the aircraft and the avionics may worsen the effect of aeromagnetic compensation. Further, according to the upper limit of the magnetic compensation, it is possible to determine whether a magnetic noise level in the area or of the aircraft is suitable for an aeromagnetic survey flight.
    [0063] It should be illustrated that the compensation coefficients in the formulas (11), (15) and (19) are real test data obtained from a certain aeromagnetic calibration flight, and then the results shown in Tables 1 and 2 are calculated.
    [0064] The compensation coefficients used in the above experiments are listed in Table 3 below.
    [0065] Table 3 Aeromagnetic compensation coefficients obtained during a calibration flight c 1 c 2 c 3 c 4 c 5 c6 1.18×10-2 1.83×10-2 4.70×10-4 1.72×10-2 -1.20×10-5 -2.55×10-5 c 7 c 8 c 9 c 10 c 11 c 12 1.68×10-2 3.19×10-5 1.62×10-2 4.5×10-3 1.1×10-3 2.76×10-4 c 13 c 14 c 15 c 16 c 17 c 18 1.64×10-4 5.1×10-3 8.28×10-4 -3.14×10-5 -5.04×10-4 4.3×10-3
    [0066] In step S106, it is determined, according to the upper limit of the error of the magnetic noise, whether the environmental background field and the aeromagnetic flight platform are suitable for an aeromagnetic survey flight, and if so, a calibration flight is performed to acquire a compensation coefficient.
    [0067] Different types of aeromagnetic experiments have different accuracy requirements. According to whether the calculated upper limit of the aeromagnetic compensation is less than a lower limit of accuracy of the aeromagnetic experiments, it is determined whether the environmental background field and the flight platform satisfy the aeromagnetic experiments. Specifically, high-accuracy aeromagnetic experiments usually require a FOM to be less than 1nT. If the calculated upper limit of the aeromagnetic compensation is greater than 1 nT, it proves that an area or an aircraft in flight is not suitable for aeromagnetic experiments under this accuracy.
    [0068] According to the above upper limit of the error, it is determined whether the environmental background field and the aeromagnetic flight platform are suitable for an aeromagnetic survey flight, and if so, a calibration flight is performed and then compensation coefficients are obtained. A least square solution of the compensation coefficients in the formula (3) is expressed as the following expression: C =A T A − 1A T H d
    [0069] The method for acquiring the compensation coefficients comprises: performing a calibration flight, and acquiring an attitude term according to outputs of the three-axis fluxgate magnetometer in three directions and strength of a geomagnetic field; acquiring an interference magnetic field of the aircraft at the optically pumped magnetometer according to data of the calibration flight; and calculating the compensation coefficients using the attitude term and the interference magnetic field.
    [0070] In the present embodiment, the attitude term in the formula (24) is obtained by using the outputs of the three-axis fluxgate magnetometer in the three directions and the strength of the geomagnetic field, and Hd is magnetic interference data which is derived through the curves obtained by the optically pumped magnetometer during the calibration flight.
    [0071] In step S108, data of an attitude term is acquired by performing an actual flight, an interference magnetic field is calculated by the data of the attitude term and the compensation coefficient acquired during the calibration flight, and magnetic compensation is performed.
    [0072] In an actual flight process, attitude information of the aircraft with respect to the geomagnetic field is measured by using the three-axis magnetometer. It may be known from the formulas (8) to (10) that the corresponding attitude term may be determined according to the direction cosines and the geomagnetic field, and therefore corresponding 18 attitude terms may be obtained by the outputs of the three-axis magnetometer and the value of the geomagnetic field, and then the interference magnetic field of the aircraft is obtained by using the 18 attitude terms and the compensation coefficients obtained during the calibration flight.
    [0073] A magnetic field after magnetic compensation is performed satisfies the following expression:H C = H M − A M C
    [0074] Since the optically pumped magnetometer has much higher sensitivity than the fluxgate, the upper limit of the aeromagnetic compensation in the formula (25) is determined by a contribution component AM C of the fluxgate.
    [0075] It should be emphasized that the magnetic compensation method based on an aeromagnetic compensation error model according to the present disclosure is not limited to the specific methods mentioned in the embodiments, and those skilled in the art may simply make changes or substitutions thereto, for example:the aeromagnetic compensation model does not necessarily comprise all 18 items as selected items to be used, and may only cover some of the 18 items.
    [0076] In summary, the present disclosure provides a magnetic compensation method based on an aeromagnetic compensation error model. Before an aeromagnetic compensation flight is performed, a model error caused by an environmental background magnetic noise, engines of an aircraft and avionics is firstly calculated to determine an upper limit of a background error, so as to determine whether a magnetic noise level in an area or of the aircraft is suitable for an aeromagnetic survey flight. There is provided a method for calculating an upper limit of magnetic compensation, which may estimate an environmental background field in an exploration area before an aeromagnetic compensation calibration flight is performed, and magnetic characteristics of an aeromagnetic flight platform are estimated effectively in advance, to determine in advance whether the exploration area and an aeromagnetic flight platform are suitable for aeromagnetic exploration, which provides a theoretical basis for the selection and construction of the platform, and avoids occurrence of an invalid flight.
    [0077] Of course, according to actual needs, the magnetic compensation method based on an aeromagnetic compensation error model according to the present disclosure further comprises other commonly-used methods and steps. Since the other commonly-used methods and steps have nothing to do with the innovation of the present disclosure, they will not be repeated here.
    [0078] The specific embodiments described above further describe the purposes, technical solutions, and beneficial effects of the present disclosure in detail. It should be understood that the above description is only specific embodiments of the present disclosure and is not intended to limit the present disclosure. Any modifications, equivalent substitutions, or improvements etc. made within the spirit and principle of the present disclosure, shall be included in the protection scope of the present disclosure.
    权利要求:
    Claims (10)
    [0001] A magnetic compensation method based on an aeromagnetic compensation error model, comprising:
    acquiring an upper limit of an error of a magnetic noise by using the aeromagnetic compensation error model, before an aeromagnetic flight, wherein the magnetic noise is caused by both an environmental background field in an exploration area and an aeromagnetic flight platform;
    determining, according to the upper limit of the error of the magnetic noise, whether the environmental background field and the aeromagnetic flight platform are suitable for an aeromagnetic survey flight, and if so, performing a calibration flight to acquire a compensation coefficient; and
    acquiring data of an attitude term by performing an actual flight, calculating an interference magnetic field by the data of the attitude term and the compensation coefficient acquired during the calibration flight, and performing magnetic compensation.
    [0002] The magnetic compensation method according to claim 1, wherein the aeromagnetic compensation error model is established by:acquiring various interference magnetic fields, and acquiring errors of magnetic interferences corresponding to the interference magnetic fields based on expressions of the respective interference magnetic fields, so as to acquire a total magnetic interference error.
    [0003] The magnetic compensation method according to claim 2, wherein the various interference magnetic fields comprise: a constant interference magnetic field, an induced interference magnetic field, and an eddy current interference magnetic field, and the total magnetic interference error in the aeromagnetic compensation error model satisfies:Error T = error P + error I + error E
    [0004] The magnetic compensation method according to claim 3, whereinthe constant interference magnetic field HPERM (t) satisfies:H PERM t = c 1 u 1 + c 2 u 2 + c 3 u 3= A 1 c 1 + A 2 c 2 + A 3 c 3 ,
    [0005] The magnetic compensation method according to claim 4, wherein the direction cosines u 1, u 2 and u 3 satisfy:u 1 = cos X t = T t ⁄ T t 2 + L t 2 + V t 2
    [0006] The magnetic compensation method according to claim 1, wherein determining, according to the upper limit of the error of the magnetic noise, whether the environmental background field and the aeromagnetic flight platform are suitable for an aeromagnetic survey flight comprises:determining whether the environmental background field and the aeromagnetic flight platform satisfy an aeromagnetic experiment, according to whether a calculated aeromagnetic compensation upper limit is less than a lower limit of accuracy for the aeromagnetic experiment, and if the aeromagnetic compensation upper limit is less than the lower limit of the accuracy for the aeromagnetic experiment, determining that it is suitable for the aeromagnetic survey flight; otherwise, determining that it is not suitable for the aeromagnetic survey flight.
    [0007] The magnetic compensation method according to claim 1, wherein performing a calibration flight to acquire a compensation coefficient comprises:
    performing the calibration flight and acquiring an attitude term according to outputs of a three-axis fluxgate magnetometer in three directions and strength of a geomagnetic field;
    acquiring an interference magnetic field of an aircraft at an optically pumped magnetometer according to data of the calibration flight; and
    calculating the compensation coefficient using the attitude term and the interference magnetic field, wherein
    the compensation coefficient satisfies: C =A T A − 1A T H d
    [0008] The magnetic compensation method according to claim 7, wherein the interference magnetic field of the aircraft at the optically pumped magnetometer satisfies:H d t = H PERM t + H IND t + H EDDY t
    [0009] The magnetic compensation method according to claim 7, wherein a magnetic field after the magnetic compensation is performed satisfies:H C = H M − A M C
    [0010] The magnetic compensation method according to any of claims 1 to 9, whereinacquiring the upper limit of the error of the magnetic noise caused by the environmental background field in the exploration area comprises: measuring the environmental background magnetic field at a certain point in the area to be explored in a state without an exploration aircraft, and acquiring the upper limit of the error of the magnetic noise caused by the environmental background field, by using the aeromagnetic compensation error model; andacquiring the upper limit of the error of the magnetic noise caused by the aeromagnetic flight platform comprises: measuring a magnetic field of the noise in an operation state, wherein, in the operation state, devices such as engines and avionics of an aircraft are turned on during an exploration operation, and acquiring the upper limit of the error of the magnetic noise caused by the aeromagnetic flight platform, by using the aeromagnetic compensation error model.
    类似技术:
    公开号 | 公开日 | 专利标题
    MacLeod et al.1993|3-D analytic signal in the interpretation of total magnetic field data at low magnetic latitudes
    CN102678106B|2015-04-01|Data processing method for LWD | electromagnetic wave resistivity logging instrument
    CN104931028B|2017-03-01|A kind of three axle magneto-electronic compass error compensation methods based on deep learning
    KR101783050B1|2017-09-28|Method and device for calibrating a three-axis magnetic field sensor
    CA2412081C|2010-01-12|System, method and computer product for geological surveying utilizing natural electromagnetic fields
    US7027923B2|2006-04-11|Method for determining sonde error for an induction or propagation tool with transverse or triaxial arrays
    US7835879B2|2010-11-16|Method and device for calibrating a magnetic sensor
    US20170176629A1|2017-06-22|Method for multi-tubular evaluation using induction measurements
    Georges et al.1970|Wave‐induced fluctuations in ionospheric electron content: A model indicating some observational biases
    Holtham et al.2010|Three-dimensional inversion of ZTEM data
    US3490032A|1970-01-13|Method and apparatus utilizing a pair of spaced magnetometers for making magnetic surveys
    AU2003202323B2|2007-11-22|Airborne vector magnetic surveys
    Hou et al.2009|A new trifilar pendulum approach to identify all inertia parameters of a rigid body or assembly
    CN105891755B|2018-07-03|The bearing calibration of aircraft hanging fluxgate magnetic gradient tensor instrument
    CN104344804B|2017-12-15|Satellite Simulation zero-g state unit pointing accuracy measuring method
    EP2820405B1|2019-05-08|Fault detection for pipelines
    CN107024674B|2019-04-26|A kind of magnetometer method for quick field calibration based on least square method of recursion
    Ferraccioli et al.1998|Microlevelling procedures applied to regional aeromagnetic data: an example from the Transantarctic Mountains |
    CN106353824B|2018-03-27|System compensation and magnetic disturbance the compensation fusion method of airborne flux-gate magnetic gradient tensor instrument
    FR2910050A1|2008-06-20|Determining the properties of a terrestrial formation using the electromagnetic coupling tenseur
    CN103837900B|2016-08-31|A kind of buried cable localization method based on Vector Magnetic Field detection and device
    AU2013227426B2|2015-10-08|Fault detection for pipelines
    US20180180417A1|2018-06-28|Systems and methods for calibrating and adjusting a heading reference system
    CN102096113A|2011-06-15|Time-domain ground-air electromagnetic detection system and calibration method
    CN106405658B|2018-03-27|A kind of campaign-styled locating magnetic objects method based on vector gradometer
    同族专利:
    公开号 | 公开日
    WO2019061403A1|2019-04-04|
    US20200301040A1|2020-09-24|
    CN111433634A|2020-07-17|
    引用文献:
    公开号 | 申请日 | 公开日 | 申请人 | 专利标题
    法律状态:
    2019-04-06| STAA| Information on the status of an ep patent application or granted ep patent|Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE |
    2020-07-03| STAA| Information on the status of an ep patent application or granted ep patent|Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE |
    2020-07-03| PUAI| Public reference made under article 153(3) epc to a published international application that has entered the european phase|Free format text: ORIGINAL CODE: 0009012 |
    2020-08-05| 17P| Request for examination filed|Effective date: 20200429 |
    2020-08-05| AK| Designated contracting states|Kind code of ref document: A1 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
    2020-08-05| AX| Request for extension of the european patent|Extension state: BA ME |
    2021-01-06| DAV| Request for validation of the european patent (deleted)|
    2021-01-06| DAX| Request for extension of the european patent (deleted)|
    2021-06-30| A4| Supplementary search report drawn up and despatched|Effective date: 20210531 |
    2021-06-30| RIC1| Information provided on ipc code assigned before grant|Ipc: G01V3/38 20060101AFI20210525BHEP Ipc: G01V3/40 20060101ALI20210525BHEP Ipc: G01V3/16 20060101ALI20210525BHEP |
    优先权:
    申请号 | 申请日 | 专利标题
    [返回顶部]